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ABSTRACT 

We analyse the abundance ratios of the light elements Mg, Ca, C and N, relative to Fe, for 147 
red-sequence galaxies in the Coma cluster and the Shapley Supercluster. The sample covers 
a six-magnitude range in luminosity, from giant ellipticals to dwarfs at M* + 4. We exploit 
the wide mass range to investigate systematic trends in the abundance ratios Mg/Fe, Ca/Fe, 
C/Fe and N/Fe. We find that each of these ratios can be well modelled using two-parameter 
relations of the form [X/Fe] = ao + a\ logo - + a2[Fe/H], where a is the velocity disper- 
sion. Analysing these "X-planes" reveals new structure in the abundance patterns, beyond the 
traditional one-parameter correlations (Mg/Fe-tr, etc). The X-planes for the a elements, Mg 
and Ca, indicate a positive correlation with velocity dispersion, and simultaneously an anti- 
correlation with Fe/H (i.e. a\ > and 02 < 0). Taking both effects into account dramatically 
reduces the scatter, compared to the traditional X/Fe— a relations. For C and N, a similar cor- 
relation with velocity dispersion is recovered, but there is no additional dependence on Fe/H 
(i.e. 01 > and 0,2 ~ 0). The explicit dependence of X/Fe on two parameters is evidence that 
at least two physical processes are at work in setting the abundance patterns. The Fe/H depen- 
dence of Mg/Fe and Ca/Fe, at fixed a, may result from different durations of star formation, 
from galaxy to galaxy, leading to a range in the degree of SNIa contribution to the chemical 
enrichment. The absence of corresponding Fe/H dependence for C and N is consistent with 
these elements being generated in intermediate and low-mass stars, on time-scales similar to 
the Fe production. The a dependence, at fixed Fe/H, is similar for elements Mg, C and N, 
despite their likely origin in stars of very different masses and lifetimes. Ca/Fe is positively 
correlated with a, at fixed Fe/H, but its dependence is significantly less steep than that of Mg, 
C and N. We discuss possible origins for this pattern of trends and find no simple prescrip- 
tion that reproduces the observations. In particular, if C and N are produced on time-scales 
comparable to that of Fe, then the X/Fe— a trends cannot be explained solely by a systematic 
variation of star-formation time-scale with a. 

Key words: galaxies: dwarf — galaxies: elliptical and lenticular, cD — galaxies: abundances 
— galaxies: evolution 



1 INTRODUCTION 

As galaxies evolve, each generation of stars releases newly synthe- 
sised elements for incorporation into later populations. The final 
pattern of elemental abundances depends on many factors such as 
the star-formation history, the stellar initial mass function (IMF), 
the yields from different masses of stars in supemovae (SNe) and 
stellar winds, the efficiency of galactic winds in removing enriched 
material from the galaxy, and the inflow of pristine or pre-enriched 
material from an external reservoir. Sophisticated chemical evolu- 
tion models have been developed using proposals for these ingredi- 
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ents, which are all uncertain to some degree (e.g. Matteucci 2001, 
Gibson et al. 2003 and references therein). 

In our Galaxy and other resolved Local Group systems, the 
models are constrained by the run of abundance ratios (e.g. X/Fe for 
element X) with metallicity Fe/tQ measured for individual stars 
from high-resolution spectra (Mc William 1997). Chemical evolu- 
tion modelling for a wide variety of formation histories is desir- 
able to test the model ingredients over the full range of conditions 



1 We use the term metallicity interchangeably with iron abundance Fe/H 
throughout the paper. We will never refer to "total" metallicity, since it is 
dominated by O, the abundance of which is poorly constrained for external 
galaxies. 
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relevant to galaxy formation. For example, if the IMF is different 
in massive starbursts than in quiescent star-forming environments 
(e.g. as proposed by Baugh et al. 2005), this should be revealed in 
different enrichment patterns. Since not all classes of galaxy are 
represented among the systems amenable to resolved stellar spec- 
troscopy, with early-type galaxies notably absent, it is desirable 
also to compare models with data for more distant galaxies, where 
only luminosity-weighted average abundances can be measured. 

Age, metallicity and abundance ratio measurements for unre- 
solved galaxies can be obtained from intermediate-resolution spec- 
tral indices, carefully calibrated via stellar libraries and isochrone 
synthesis models (e.g. Worthey et al. 1994; Vazdekis et al. 1999; 
Thomas, Maraston & Bender 2003a; Schiavon 2007). A key mile- 
stone in developing such models was the inclusion of non-solar 
abundances of Mg (or more generally the a elements), first by Tan- 
talo, Chiosi & Bressan (1998) and Trager et al. (2000a) and then in 
the widely-used models of Thomas and collaborators (Thomas et 
al. 2003a; Thomas, Maraston & Korn 2004). The same basic ma- 
chinery has since been generalized to allow for more flexible ele- 
mental mixtures (e.g. Thomas, Maraston & Bender 2003b; Schi- 
avon 2007). Among the other light elements which are readily 
measurable in unresolved spectra are Ca, C and N (e.g. Graves & 
Schiavon 2008). Serven, Worthey & Briley (2005) argue that many 
more elements could be measured from high signal-to-noise spec- 
tra, using an enlarged set of index definitions. 

These analysis techniques were initially applied to samples 
comprised mainly of giant elliptical and SO galaxies (e.g. Trager et 
al. 2000b; Kuntschner et al. 2001 ; Thomas et al. 2005). Later work 
extended the coverage to a wider range of luminosity (e.g. Cald- 
well et al. 2003; Nelan et al. 2005; Smith, Lucey & Hudson 2007). 
The increased baseline helped to reveal the overall scaling relations 
of age, metallicity and Mg/Fe, as a function of velocity dispersion, 
a. (Stellar population parameters are generally more tightly related 
to a than to other mass proxies such as luminosity.) On average, 
less massive galaxies are younger, less metal-rich and have lower 
Mg/Fe, than the giant ellipticals. The Mg/Fe— a trend has been in- 
terpreted as evidence for a shorter time-scale of star formation in 
high-mass galaxies, preventing the incorporation of Fe-rich ejecta 
from Type la supernovae (SNIa), which are delayed relative to star 
formation, into the stars observed today (Thomas et al. 2005). Other 
scenarios have been discussed, e.g. by Trager et al. (2000b), includ- 
ing an IMF which varies as a function of galaxy mass, altering the 
relative number of SNIa and SNII, or selective galactic winds to 
modulate the effective yields. 

For other light elements, the picture is less clear. Ca was ex- 
pected to be enhanced similarly to Mg, since both are produced in 
the triple-Q chain and released in Type II supernovae (SNII). In- 
stead, Ca appears to have only solar abundance with respect to Fe 
in giant ellipticals (Cenarro et al. 2003; Thomas et al. 2003b). For 
the abundances of C and N, which are probably produced mainly 
by less massive stars, somewhat contradictory results have been ob- 
tained (e.g. Clemens et al. 2006; Kelson et al. 2006; Graves et al. 
2007). 

In this paper, we analyse abundance ratios measured for 81 
dwarf galaxies in the Coma cluster (from Smith et al. 2009), and 
66 generally more luminous galaxies in the Shapley Supercluster 
(Smith et al. 2007). The key observational improvement over pre- 
vious work is that abundance ratios are determined for individual 
galaxies, rather than ensemble averages, and we probe a very wide 
range in luminosity, velocity dispersion and metallicity. A second 
advance is the introduction of planar relationships between the 
abundance ratios, the metallicity and the velocity dispersion. These 



new scaling relations are potentially powerful constraints on the 
ingredients of chemical evolution models for unresolved galaxies. 
The outline of the paper is as follows. Section|2]briefly reviews the 
properties of the two galaxy samples and the data employed, and 
describes the fits we use to model the abundance ratio correlations. 
In Section [3] we investigate empirically the correlations between 
the abundance ratios and the luminosity, velocity dispersion, age 
and metallicity of the galaxies. We discuss in turn the behaviour of 
Mg, Ca, C and N, in the context of previous observational work. 
Section[4]combines the information obtained for the four elements, 
and explores how the results relate to models for chemical evolu- 
tion. Our principal conclusions are summarized in Section|5] 



2 DATA AND METHODS 

The data used in this paper are based on spectroscopy from Hec- 
tospec at the 6.5m MMT for dwarf galaxies in Coma (Smith et 
al. 2009), and from AAOmega at the 3.9m Anglo-Australian Tele- 
scope for more luminous galaxies in the Shapley Supercluster 
(Smith et al. 2007). The absorption line index data are tabulated 
in these papers. Both studies targeted galaxies in the red sequence, 
through explicit colour selection (Coma), or post-facto removal of 
galaxies with Ha emission (Shapley). The treatment of the two 
datasets is very similar in all respects. Velocity dispersion mea- 
surements are available for 57 of 66 Shapley galaxies, from the 
AAOmega observations, and for 42 of 81 Coma galaxies, compiled 
from literature sources. The line indices were used to derive stellar 
population parameters by comparison to the models of Schiavon 
(2007), using the EZ-AGES code (Graves & Schiavon 2008). The 
parameters obtained are SSP-equivalent age, metallicity Fe/H and 
the light-element abundance ratios Mg/Fe, Ca/Fe, C/Fe and N/Fe. 
The Mg abundance is derived from the Mgb5177 index, Ca from 
Ca4227, C from C4668 and N from the CN band at 4175 A. Esti- 
mates of the parameter errors and covariances were derived from a 
Monte-Carlo error model. Full details of this process, and tables of 
parameter estimates, are presented by Smith et al. (2009). All quan- 
tities used in this paper refer to measurements within an aperture of 
physical radius 0.95 kpc (~0.5-1.0x the effective radius, i? c ff) in 
Shapley, and 0.35 kpc (~0.2-0.4xi? cff ) in Coma. 

Passive galaxies follow many inter-related correlations, such 
as the so-called "Z-plane" (Trager et al. 2000b), which relates ve- 
locity dispersion, SSP-equivalent age and metallicity Fe/H. A fit to 
the Z-plane for our combined Coma and Shapley sample yields 

[Fe/H] = (0.69 ± 0.07) log a - (0.70 ± 0.09) logt SS p - 0.90 , 

where the coefficients are all within la of the equivalent plane from 
Trager et al. Although abundance ratios such as Mg/Fe have tradi- 
tionally been modelled as functions of a, it is not clear a priori 
whether the fundamental dependence is on velocity dispersion, or 
on age, metallicity or some other related parameter. The wide range 
in galaxy properties spanned by our combined sample will be ex- 
ploited in this paper to disentangle the effects on X/Fe of these 
mutually-correlated properties. 

We will fit models simultaneously to the two samples, in or- 
der to cover the widest baseline in mass and metallicity. The Coma 
sample provides excellent coverage of the dwarf galaxy regime, 
while the Shapley data improve the leverage at high a. Qualitatively 
similar results, albeit with larger errors, can be derived from either 
sample taken individually (see Appendix [A). As a first step, we fit 
simple correlations of X/Fe with luminosity Lr, velocity disper- 
sion a, SSP-equivalent stellar age issp and metallicity Fe/H, con- 
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Figure 1. Correlations of Mg/Fe with luminosity, velocity dispersion (where available), age and metallicity Fe/H. The Coma galaxies are shown with blue 
filled symbols and the Shapley comparison galaxies with red open symbols. Indivdual error bars are shown for the luminosity and velocity dispersion panels 
only; in the age and metallicity panels, we show a typical error ellipse for each sample in the lower left. Lines show fits to the combined sample (dashed black), 
to Coma only (bold blue) and to Shapley only (light red). The yellow track in the velocity dispersion panel indicates the mean trend for giant galaxies, from 
Graves et al. (2007). The title-bar gives the measured slope, and total rms scatter for the fit to the combined sample (see Table [T] for fits to the Coma and 
Shapley data separately, see Table lAlV The product-moment correlation coefficient is r. The probability, /, of exceeding this r-value by chance is computed 
from simulations which include the effect of correlated errors (see Section[2]for details). The r expected from correlated errors alone is given in parentheses. 
Age and velocity dispersion are the best single-parameter predictors of the Mg/Fe ratio; a simultaneous fit prefers a combination of a and Fe/H dependences. 
Note that the "offset" between Shapley and Coma in the Fe/H panel arises because the Shapley galaxies have much higher a on average. 




log a log a log a 



Figure 2. Left: the Fe/H-cr correlation for the combined Shapley (red) and Coma (blue) sample. Although these parameters are correlated, our sample covers 
sufficient baseline in each parameter that correlations primarily with a can be distinguished from correlations primarily with Fe/H. Centre: the traditional 
Mg/Fe-cr relation. Right: the cr-projection of the X-plane for Mg. The vertical errorbars account for the anti-correlated errors in Mg/Fe and Fe/H. Introducing 
the additional Fe/H dependence reduces the total rms scatter is from 0. 10 dex to 0.07 dex (a factor of two in the variance), and the intrinsic scatter from 0.07 dex 
to 0.04 dex (a factor of three in variance). Note that the reduction of scatter can be appreciated in both the Shapley and the Coma sample taken individually. 
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Figure 3. Correlations for Ca/Fe. The behaviour is qualitatively similar to that of Mg/Fe in Figure [T] except for the absence of correlation with a. Symbols 
and annotations as in Figure [T] 



sidering each of the latter quantities in turn. These one-parameter 
tests are summarized in Table Q] 

For the correlations with age and metallicity, which are esti- 
mated simultaneously with the abundance ratios in the EZ-AGES 
process, it is important to consider the effect of the correlated er- 
rors. To this end we compare the product-moment correlation co- 
efficient r of the observed data against 100 synthetic datasets for 
each panel. The simulations assume the "true" Xi (e.g. [Fe/H]) are 
as observed, while the true Yi (e.g. [Mg/Fe]) are zero. The fake "ob- 
served" data (Xi,yi) are drawn from the error ellipse, and the corre- 
lation coefficient recorded. The fraction / of simulations yielding 
r greater than the observed (one-tailed) is then a measure of the 
significance which accounts for the correlated errors. In fact, the 
errors are sufficiently small that the distribution of galaxies in Fe/H 
(and age) is well-resolved. For example, the observed rms in Fe/H 
values is 0.24 dex, while the errors can account for only O.lOdex. 
As a result, in all cases where apparently significant trends are re- 
covered, the simulation tests confirm that they cannot be generated 
spuriously by the correlation of errors. 

Considering the one-parameter fits, especially those for Mg 
(Section |3.1t , it became clear that introducing a second parame- 
ter would provide a significantly improved description of the X/Fe 
variations. The two-parameter model with smallest scatter uses the 
velocity dispersion a and metallicity Fe/H as predictors of the 
X/Fe, in the form: 

[X/Fe] = a + a, log ( — £ - ) + a 2 [Fe/H] . 
VlOOkms / 

These fits are restricted to the 99 galaxies in the combined sample 



which have velocity dispersion measurements. We prefer to em- 
ploy a rather than luminosity here, because it is not affected by 
variations in the stellar mass-to-light ratio, and because in practice 
tighter correlations are obtained. The coefficients and scatters of 
these relationships, which we call the "X-planes", are summarized 
in Table|2] Although we have not explicitly accounted for the corre- 
lated measurement errors in fitting the planes, we do allow for this 
effect when estimating the intrinsic scatter around the X-planes. 
Moreover, the simulations used to derive the scatter estimates also 
yield estimates for the coefficient bias due to the error correlations, 
which is confirmed to be minimal (~0.01 for Mg and Ca, ~0.04 
for C and N). 

For completeness, Table [3} presents the results of more gen- 
eral fits, of the form 



[X/Fc] = b + &i log 



100 kms~ 



- + b 2 [Fe/H] + 



63 log 



^SSP 

lOGyr 



+ b 4 S . 



where the fourth term allows for additional variations with SSP- 
equivalent age, tssp- In the fifth term, S = 1 for Shapley galaxies 
and S = for Coma galaxies. Its coefficient b 4 is treated as a nui- 
sance parameter to account for unknown offsets between the two 
data sources, including aperture effects. We do not include lumi- 
nosity as a further predictor here, since it is highly redundant with 
a. (Experiments with including luminosity do not suggest any sig- 
nificant extra dependences.) In general, we find that the coefficients 
&3 and 64 above are consistent with zero, and that using the more 
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Figure 4. Correlations for Ca/Mg. This figure highlights the differences in behaviour for these two a elements, in particular the anti-correlation with velocity 
dispersion. Symbols and annotations as in Figure[T] 



general fits does not reduce the residual scatter significantly. This 
justifies, post facto, the use of the two-parameter X-planes as a con- 
cise description of the abundance correlations. 

Note that the X-plane incorporates velocity dispersion and 
metallicity on an equal footing as statistical predictors of the abun- 
dance ratios. However, while a can plausibly be interpreted as a 
causal predictor of the abundances (i.e. it is an input to chemical 
evolution models), this is not true for Fe/H, which should be an 
output of any chemical evolution model. Instead, a correlation of 
X/Fe with Fe/H must indicate a mutual dependence of Fe/H and 
X/Fe on some other aspect of the formation process. We discuss 
this further in Section|4] 



3 ABUNDANCE RATIO TRENDS 

In this section, we investigate empirically the correlations followed 
by the abundance ratios Mg/Fe, C/Fe, N/Fe and Ca/Fe. We will also 
make use of some ratios between the light elements, specifically 
Ca/Mg and C/N, to isolate similarities and differences among pairs 
of elements with similar expected origin. 

3.1 Magnesium 

The overabundance of Mg relative to Fe in giant elliptical galax- 
ies has long been recognized (e.g. Peterson 1976). The systematic 
trends in Mg/Fe (or a/Fe, nearly equivalently) have been investi- 
gated many times, especially since the publication of the Thomas 
et al. (2003a) models (Trager et al. 2000b; Kuntschner et al. 2001; 



Thomas et al. 2005; Nelan et al. 2005 and many others). A consen- 
sus has emerged that Mg/Fe rises with velocity dispersiorQ reach- 
ing ~0.3 dex above solar for giant ellipticals with a >200 km s _1 . 
For dwarf galaxies, the reported Mg/Fe ratios are typically solar 
(Geha et al. 2003) or modestly super-solar (e.g. Chilingarian et al. 
2008). A hint that more complex correlations may be present was 
found by Sanchez-Blazquez et al. (2006), who showed that at fixed 
a, older galaxies have higher Mg/Fe on average (their Figure 11). 

The single-parameter correlations for Mg/Fe in our sample are 
presented in Figure Q] We recover the usual increase with velocity 
dispersion, and a weaker correlation with luminosity. The Coma 
and Shapley samples seem to follow a common Mg/Fe-a rela- 
tion, the slope of which is similar to that of recent work, including 
Graves et al. (2007). At the lowest a covered, the average Mg/Fe re- 
mains slightly enhanced over solar (as found by Chilingarian et al. 
2008). There is also a correlation with age, £ssp, which is in fact 
stronger than the a trend. Again, the Coma and Shapley galaxies 
follow a similar relation. 

A more complex correlation structure is revealed in the panel 
showing Mg/Fe as a function of Fe/H. Here we might naively 
have expected a positive trend, since high-Fe/H galaxies tend to 
be more massive, and more massive galaxies have enhanced Mg/Fe 



2 A superficially contrary result was obtained by Kelson et al. (2006), who 
found a/Fe flat or even decreasing with a in a distant galaxy cluster. Note 
however, that the usual Mgb5177 index was not used in their fits, hence it is 
not clear which a elements were being measured, and a direct comparison 
with Mg/Fe may not be justified. 
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Figure 5. Correlations for C/Fe. Note the steep and tight correlation with velocity dispersion, and the absence of dependence on Fe/H. The latter is quite unlike 
the case for the a elements (Figures[T]and|3j. Symbols and annotations as in Figure[T] 



(Figure [2] left and centre). Instead, we recover a significant anti- 
correlation of Mg/Fe with Fe/H. (The Monte Carlo simulations 
confirm that this result cannot be generated from the correlated er- 
rors, at the > 95% confidence level.) The Coma and Shapley sam- 
ples follow parallel but offset relations. Since the Shapley galax- 
ies have higher velocity dispersion on average, this suggests that 
Mg/Fe depends simultaneously on both a and Fe/H, and that a 
single-parameter model is insufficient to capture this complex be- 
haviour. 



The simultaneous fits (Tables[2]and[3j confirm that Mg/Fe de- 
pends on a linear combination of Fe/H and velocity dispersion. On 
average, a galaxy will have a large over-abundance of Mg if its ve- 
locity dispersion is large, but also if its metallicity is low. In the 
X-plane, the slope of the a dependence is much steeper than in 
the single-parameter Mg/Fe— a relation. This is because the high-cr 
galaxies tend to have high Fe/H on average, which partially com- 
pensates for the higher velocity dispersion. The X-plane provides a 
factor-of-two reduction in variance, compared to the Mg/Fe-a re- 
lation (Figure[2] centre and right). The observed rms scatter around 
the X-plane is only ~0.07 dex, most of which can be accounted for 
by the estimated errors (allowing for correlated errors). The esti- 
mated intrinsic component is only ~0.04 dex. The fairly tight (one- 
parameter) Mg/Fe-tssp relation seems to be a consequence of the 
X-plane correlation, coupled with the anti-correlation of age and 
Fe/H at fixed a (the Z-plane of Trager et al. 2000b). The X-plane is 
much better predictor of Mg/Fe than the Mg/Fe-tssp relation, and 
no additional age correlation is recovered in the full simultaneous 



fit. The X-plane is also much tighter than a two-parameter model 
based on a and fssp- 

Note that although our default results are based on the com- 
bined Shapley and Coma sample, fitting the two samples separately 
yields consistent coefficients for the X-planes. From each sample 
we recover significant positive correlation with a and significant 
anti-correlation with Fe/H. This confirms that the X-plane is not 
driven by "offsets" between the two samples, but is a continuous 
correlation followed consistently by galaxies in both clusters. The 
results of these separate-sample fits are provided in Appendix [Al 

We conclude that the X-plane provides the most accurate two- 
parameter description of the Mg/Fe data, is much preferred over 
one-parameter models, and is not significantly improved upon by 
fits with more parameters. Motivated by this finding, we continue 
by applying the same form of model to the other measured abun- 
dance ratios. 



3.2 Calcium 

It has long been suggested (e.g. Vazdekis et al. 1997; Worthey 
1998) that the abundance of Ca does not behave identically to that 
of Mg, despite the expectation that the production of both of these 
a-elements should be dominated by SN II. Thomas, Maraston & 
Bender (2003b) showed that Ca/Mg is subsolar by ~0.2dex for 
giant ellipticals with <r>250kms~ 1 , but that Ca/Mg is approxi- 
mately solar for low-mass galaxies a < 100 km s - , including Lo- 
cal Group dwarf spheroidals. In the Lick system, Ca is traced by the 
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Figure 6. Correlations for N/Fe. The trends are strikingly similar to those for C/Fe (Figure[5j. Symbols and annotations as in Figure[T] 



Ca4227 indejfl which is sensitive to the neighbouring CN band- 
head, as discussed by Prochaska, Rose & Schiavon (2005). The 
EZ-AGES approach in principle allows for the effect of varying 
C/Fe and N/Fe, and should help to decouple the Ca abundances. 
Even with this method, however, Graves et al. (2007) report only a 
weak increase of Ca/Fe with a for giant red-sequence galaxies, with 
approximately solar abundance on average, suggesting the anoma- 
lous results do not arise from CN contamination. 

Our correlations for Ca/Fe are shown in Figure[3] We confirm 
the absence of a significant trend with velocity dispersion, or with 
luminosity, although our average Ca/Fe ratios are slightly super- 
solar throughout. There is however a strong anti-correlation with 
Fe/H, which is very similar to the case for Mg/Fe, and a weak cor- 
relation with age is also recovered. 

The X-plane and multi-parameter fits confirm that the domi- 
nant dependence for Ca/Fe is an anti-correlation with Fe/H. With 
the simultaneous fits, we do recover an increase of Ca/Fe with <j, 
at fixed Fe/H. This was "hidden" in the single-parameter fits by the 
metallicity trend: high-cr galaxies tend also to have high Fe/H, caus- 
ing compensating effects on Ca/Fe. In the multi-parameter fit, there 
is marginal evidence for additional correlation of Ca/Fe with age, 
although this makes insignificant contribution to the total variance. 

To highlight the similarities and differences between Mg and 
Ca, we show in Figure [4] the equivalent correlations for Ca/Mg. 
The strongest trend is an anti-correlation with velocity dispersion. 



3 An alternative tracer is provided by the Ca triplet at 8500-8700 A, which 
has also been used to infer under-abundant Ca, and/or variations in the stel- 
lar initial mass function (Saglia et al. 2002; Cenarro et al. 2003). 



The Shapley and Coma samples seem to follow a common relation, 
with Ca under-abundant by ~0. 1 dex relative to Mg in galaxies with 
cr>50kms" 1 or L r > 1O 9 L , but near solar in dwarf galaxies. 
We conclude that Ca indeed behaves somewhat differently from Mg 
in the X-plane formulation. However, the relations for Ca and Mg 
share similar features when viewed in contrast to the very different 
correlations found for elements C and N which we explore next. 



3.3 Carbon 

There is a long-standing debate as to whether short-lived high-mass 
stars (via Wolf-Rayet winds) or long-lifetime intermediate/low- 
mass stars (through the planetary nebula phase) are primarily re- 
sponsible for enrichment of C in the solar neighbourhood (e.g. 
Gustafsson et al. 1999; Carigi et al. 2005). 

Determinations of C/Fe, and its trend with mass, in unresolved 
galaxies have been somewhat contradictory. Sanchez-Blazquez et 
al. (2006) report C/Fe rising with a, but less steeply than Mg, while 
Clemens et al. (2006) find C/H (and therefore also C/Fe) declining 
with a, and Kelson et al. (2006) argue for no strong a dependence 
in C/Fe. More recently, Graves et al. (2007), using EZ-AGES, re- 
cover a C/Fe increase with a that is even steeper than that of Mg/Fe. 

Figure[5] shows the correlations followed by C/Fe in our data. 
We find that C/Fe increases with luminosity and with velocity dis- 
persion. The Shapley and Coma galaxies follow a similar corre- 
lation with a, where the slope is also consistent that reported by 
Graves et al. for giant galaxies. In the luminosity correlation, by 
contrast, there is some evidence for a change in slope between gi- 
ants and dwarfs. There is a weak increase of C/Fe with age (as seen 
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Figure 7. Correlations for C/N. The absence of significant trends in this figure emphasizes the systematic similarity in the behaviour of C and N. Symbols and 
annotations as in Figure[T] 



for both Mg/Fe and Ca/Fe). However, the behaviour of C differs 
strikingly from that of Mg and Ca in the metallicity panel: we do 
not recover a strong anti-correlation of C/Fe with Fe/H, and indeed 
the correlation is marginally positive. 

The simultaneous fits again confirm the above discussion. The 
behaviour of C/Fe can be described by a strong underlying correla- 
tion with a. After accounting for this, there is a marginally signifi- 
cant anti-correlation with Fe/H, which was hidden by the a depen- 
dence. However, the metallicity trend of C is very much shallower 
than the equivalent trends for Mg and Ca. 



Graves et al. (2007) report an increase in N that is much stronger 
than the equivalent for MjjJ 

Our results for N are presented in Figure [6] and are strikingly 
similar to the equivalent relations for C: there are strong correla- 
tions of N/Fe with luminosity and velocity dispersion, and only 
weak dependence on age and metallicity. Our N/Fe— a slope is sim- 
ilar to that reported by Graves et al. for giant galaxies, but our rela- 
tion has a higher zero-point, by ~0. 1 dex. 

As expected from Figure [6] the coefficients of the X-plane for 
N are indistinguishable from those of C. To emphasise the absence 
of systematic variation between these elements, Figure [7] shows 
single-parameter correlations for the C/N ratio, where no signifi- 
cant slopes are obtained. 



3.4 Nitrogen 

The production of N is thought to be dominated by intermediate- 
mass stars, through mass loss in the asymptotic giant branch (AGB) 
and planetary nebula phases (van den Hoek & Groenewegen 1997), 
although there could also be contributions from rapidly-rotating 
massive stars (Meynet & Maeder 2002). Previous work on early- 
type galaxies has suggested that N/Fe is enhanced over solar, and 
that the enhancement increases with a. Sanchez-Blazquez et al. 
(2006) found that to reproduce their observed index-a relations, 
N/Fe should increase with a but less steeply than Mg/Fe. Kelson et 
al. (2006) found N/Fe rising steeply with velocity dispersion (con- 
verted from their quoted a/N and afFe trends). Using EZ-AGES, 



4 DISCUSSION 

As noted in the previous section, our results for the average trends 
of X/Fe as a function of velocity dispersion largely agree with 
previous work in the regime of giant galaxies. In particular, our 
(single-parameter) X/Fe— a slopes are consistent with those of 
Graves et al. (2007), who used the the same stellar population 
models (Schiavon 2007) and inversion method (Graves & Schiavon 
2008). 

4 Very recently, Toloba et al. 2008 report only a very weak dependence of 
the NH 3360 A band on velocity dispersion, which appears contrary to the 
results from the CN band used here and in most other work. 
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Figure 8. Confidence intervals on the parameters of the X/Fe planes. For 
each element, the ellipses show the 68% and 95% confidence intervals on 
the two slope coefficients. The dashed lines indicate zero variation, while 
the dotted lines coincide with the parameters for the Mg plane. 

Where our study breaks new ground is in probing a much 
larger range in galaxy mass and metallicity, and in analysing the 
abundance patterns using multi-parameter fits. The results of our 
X-plane fits for all four elements are summarized in Figure [8] in 
the form of confidence intervals on the fit coefficients. (An equiv- 
alent figure for the separate X-plane fits to Coma and Shapley is 
shown in Appendix [A] and confirms that the main results can be 
obtained from either sample individually.) An alternative way to 
visualise the bivariate dependence of X/Fe is by means of contours 
in the Fe/H-a plane, as shown in Figure [9] The X-planes are read- 
ily seen in our data owing to the wide baseline covered in both a 
and Fe/H. These results could not easily have been recovered from 
earlier samples which were dominated by metal-rich high-a giant 
ellipticals. 

A detailed comparison against realistic chemical evolution 
models is clearly required to interpret the results we have presented. 
In the following sections, we limit the discussion to some qual- 
itative arguments. First, we argue that variation in star-formation 
time-scale, at fixed a, is responsible for imprinting the a/Fe-vs- 
Fe/H trends, while leaving the C/Fe and N/Fe ratios independent of 
Fe/H (Section [4.1t . Second the comparable increase of X/Fe with 
a for all four light elements requires a separate explanation, since 
time-scale and other progenitor-mass effects should not distinguish 
between SNIa products (Fe) versus AGB products (C and N) (Sec- 
tion [472} . Third, we comment on the apparent discrepancy in the 
behaviour of Ca compared to Mg (Section[43}- 

4.1 Fe/H trends: a range in star-formation time-scale at fixed 
velocity dispersion 

For individual stars in the solar neighbourhood, the Mg/Fe ratio is 
enhanced at low Fe/H, where SNII enrichment is important, and 
drifts towards solar values at high Fe/H, as SNIa become dom- 
inant (McWilliam 1997 and references therein). This pattern re- 
sults from SNIa being on average significantly delayed relative 



to star formation (and SNII), as expected if their progenitors are 
low/intermediate mass starjf] 

Of course, in galaxies we observe only an average over the en- 
richment history, rather than the "snapshots" of chemical evolution 
revealed by individual stars. None the less, if other galaxies expe- 
rienced similar enrichment processes, operating over a variety of 
time-scales, we might expect the integrated abundances to follow a 
qualitatively similar trend. Galaxies where star formation persisted 
for longer would have more time to incorporate SNIa ejecta, raising 
the average Fe/H and lowering the average Mg/Fe ratio. 

A key puzzle in earlier work on elliptical galaxies was that 
higher-metallicity galaxies have higher Mg/Fe ratios, i.e. opposite 
to the pattern in galactic stars. By analysing the X-plane for Mg, 
we have shown that there is an anti-correlation like that in the stars, 
after referencing to a common velocity dispersion. Intriguingly, we 
find that the slope of the Mg/Fe-vs-Fe/H dependence at fixed o 
(—0.33) agrees quantitatively with that recovered for galactic disk 
stars, over a similar metallicity range (Figure [10] top panel). The 
resemblance is particularly striking for the thick disc stars taken 
alone. We may interpret this slope as the characteristic signature 
of star-formation being quenched during the transition from SNII- 
dominated to SNIa-dominated abundance patterns. (The thin disk 
under went the same transition at lower metallicity). 

For Ca, we found an X-plane which to first order resembles 
that of Mg, as expected if they are both produced mainly in SNII, 
and in particular, the anti-correlation with Fe/H is similar to that for 
Mg. Thus both the a elements qualitatively support a picture where, 
at a given velocity dispersion, different galaxies experienced differ- 
ent time-scales of star formation, leading to an anti-correlation of 
a/Fe with Fe/H. Again, the Fe/H coefficient of the X-plane is well 
matched to the Ca/Fe-Fe/H relation for thick disk stars (FigurefTOl 
second panel). 

While it is generally agreed that the a elements Mg and Ca 
are released in SNII from high-mass stars (> 8 Mq), there is less 
consensus on the stars which dominate production of C and N. In 
the case of C, there is continuing debate as to the relative contribu- 
tions from high-mass (>8 Mq) vs low/intermediate mass stars (1- 
8 M Q ) (Gustafsson et al. 1999; Carigi et al. 2005). For N, there is 
an additional dimension to the debate, on the relative importance of 
primary and secondary synthesis, i.e. whether yields depend on the 
metallicity of the gas from which the star was born (e.g. Gavilan, 
Molla & Buell 2006). 

Recent work on galactic disk stars has shown little if any in- 
crease in C/Fe and N/Fe with declining Fe/H (Israelian et al. 2004; 
Bensby & Feltzing 2006). The simplest interpretation is that the 
production of C and N keeps pace with that of Fe, i.e. these ele- 
ments are produced in lower-mass stars with similar characteristic 
lifetimes as the SNIa progenitors. Our X-plane results suggest that 
the production sites of C and N in other galaxies are likewise dom- 
inated by low/intermediate-mass stars, which would lead naturally 
to a corresponding independence of C/Fe and N/Fe with respect 
to Fe/H. Figure [TO] confirms the similarity of the disk abundance 
pattern to the X-plane slopes for C. For N, the stellar data (from a 
different literature source than the other elements) suggest a mildly 
positive correlation which is not seen in the X-plane slope. 



5 Recent evidence that a fraction of SNIa occur promptly after star forma- 
tion (e.g. Sullivan et al. 2006) does not affect this general scenario. 



© 2007 RAS, MNRAS 000,[THT3] 



10 Russell J. Smith et al. 



4.2 Trends of X/Fe with velocity dispersion 

Having discussed the variation in X/Fe at fixed a, we now con- 
sider the systematic trends as a function of a. The key constraint is 
that the a dependences of Mg, C and N are consistent with being 
identical, with X/Fe oca 0,4 at constant Fe/H. (Note that the similar- 
ity in the a trends is apparent in the X-planes, where the different 
dependences on Fe/H are accounted for, but not in the traditional 
single-parameter scaling relations.) 

Several different mechanisms have been proposed to explain 
the traditional Mg/Fe— a correlation, by altering the relative frac- 
tions of SNIa vs SNII enrichment (see discussion by Trager et 
al. 2000b). The most widespread interpretation is that galaxies of 
higher velocity dispersion formed their stars on more rapid time- 
scales (Thomas et al. 2005), so that SNIa make less contribution to 
enriching the stellar population. However, if C and N are produced 
on nearly the same time-scale as Fe, then this scenario predicts that 
only the a elements Mg and Ca should be strongly correlated with 
a. Thus the steep C/Fe— a and N/Fe— a relations (in combination 
with the weak metallicity dependence of these elements) argues 
against a systematic variation in star-formation time-scale with ve- 
locity dispersion. 

An alternative explanation for the Mg/Fe— <r relation invokes 
variation of the IMF slope in favour of higher-mass stars at high 
velocity dispersion. This would produce a higher ratio of SNII en- 
richment (Mg, Ca) compared to SNIa (Fe), as required, but would 
not increase the C and N abundances, if these element are indeed 
generated in stars of mass comparable to that of SNIa progenitors. 
Thus the steep C/Fe— a and N/Fe— a relations also argue against 
an explanation in terms of IMF variation. Selective outflows of gas 
provide another way to vary the effective yields of SNII and SNIa. 
However, any mechanism that operates only in the early SNII- 
dominated phase will fail to match the observed trends because it 
will not distinguish between Fe production in SNIa and C and N 
production in AGB stars. For later galactic winds, it is plausible 
that SNIa products are ejected by the explosion energy, while AGB 
products are retained. However, this would lead to enhanced C/Fe 
and N/Fe in low-a galaxies, opposite to the sense of the observed 
trend. A final possibility is that the number of SNIa explosions (per 
unit stellar mass) may be smaller in high-er galaxies than in low-er 
objects. This could occur if larger star-star or cloud-cloud relative 
velocities in a star-forming region somehow suppress the forma- 
tion of binary SNIa progenitors. However, this model would also 
predict low Fe abundances in the hot interstellar gas, which could 
be in conflict with recent X-ray analyses (Humphrey & Buote 2006; 
Tawara et al. 2008). 

In this discussion we have focussed on the X/Fe-a relations, 
but the Fe/H-a relation itself is also relevant here. In particular, 
models in which the Mg/Fe— a relation arises from reduced incor- 
poration of Fe-rich SNIa ejecta at high a (e.g. varying time-scales, 
SNIa rate modulation) would at face value predict lower Fe/H in 
higher-cr galaxies, in contradition to the observed positive correla- 
tion (Figure[2] left). Taken together, the variation of X/Fe and Fe/H 
with a cannot be reproduced with any of the traditional models for 
explaining abundance ratio trends, at least in their simplest form. It 
is conceivable that combinations of these scenarios, e.g. SNIa sup- 
pression with an additional non-selective galactic wind, reducing 
all metals equally in low-u galaxies, might be compatible with the 
observed correlations. Quantitative chemical evolution modelling 
is clearly required to test more complex models of this kind. 



Table 1. The single-parameter fits for X/Fe. The quoted rms is the total 
observed scatter, i.e. uncorrected for measurement error, in dex. The table 
is organised for ready comparison between the four elements. 



Ordinate 



X 



zero-point 



slope 



log I 



Mg 
Ca 

C 

N 



+0.17 ±0.01 
+0.10 ± 0.02 
+0.07 + 0.01 
+0.11 + 0.01 



+0.18 + 0.04 0.10 

-0.05 + 0.05 0.13 

+0.38 ± 0.04 0.09 

+0.37 + 0.04 0.11 



log 



10 y L o 



Mg 

Ca 
C 
N 



-0.11 + 0.01 

-0.12 ± 0.01 
-0.07 + 0.01 
-0.01 + 0.01 



-0.05 + 0.02 0.11 



-0.03 ± 0.02 
-0.13 ± 0.02 
-0.12 + 0.02 



0.14 
0.11 
0.13 



loe ( * ssp 

1Uf = I 10 Gyr ) 


Mg 
Ca 
C 
N 


+0.18 ± 0.01 
+0.14 + 0.02 
+0.03 ± 0.01 
+0.07 + 0.02 


+0.28 ± 0.04 
+0.13 ± 0.05 
+0.21 ± 0.05 
+0.19 ± 0.05 


0.10 
0.14 
0.13 
0.14 


[Fe/H] 


Mg 


+0.09 ± 0.01 


-0.16 ± 0.04 


0.11 




Ca 


+0.03 ± 0.01 


-0.29 ± 0.04 


0.12 




C 


+0.02 ± 0.02 


+0.11 ± 0.04 


0.13 




N 


+0.07 + 0.02 


+0.11 ± 0.05 


0.15 



4.3 The calcium puzzle 

Finally, we address the apparently anomalous behaviour of Ca/Fe 
with respect to velocity dispersion. The coefficient of a in the X- 
plane for Ca is much smaller than the ~0.4 shared by the other three 
elements. While Ca/Mg is approximately solar at low a, the shal- 
lower trend yields surprisingly low Ca/Mg for massive galaxies, as 
has been noted before (e.g. Thomas et al. 2003b). 

Pipino & Matteucci (2004) point out that SNIa make a non- 
negligible contribution to Ca enrichment, according to the yields 
they adopt (from Nomoto et al. 1997). This results in reduced 
separation between the Ca/Fe ratios between SNII-dominated and 
SNIa-dominated enrichment. Since the ratios are quoted relative to 
the solar value, this leads to lower than expected Ca/Fe (i.e. sub- 
solar Ca/Mg), in the SNII-enriched populations which predomi- 
nate at high a. The chemical evolution models of Pipino & Mat- 
teucci (2004) predict [Ca/Mg] ~ -0.1 5 for massive galaxies, in 
good agreement with the high-cr offset in our Figure [4] Thus sce- 
narios which modulate the relative importance of SNIa and SNII as 
a function of a should, qualitatively at least, produce a shallower 
dependence for Ca/Fe than for Mg/Fe, as observed. However, the 
same argument would also predict a shallower dependence of Ca/Fe 
on Fe/H, which is not seen in our data. 

An additional anomaly is that alone among the abundance ra- 
tios considered here, Ca/Fe exhibits a significant dependence on 
SSP-equivalent age in the extended multiparametric fit (Table [3}. 
Specifically, older galaxies have lower Ca/Fe than younger galax- 
ies with the same a and Fe/H. Table [3] shows that correcting to 
common age yields a a dependence somewhat closer to the value 
for Mg (0.22+0.06 versus 0.31+0.04). We do not speculate on the 
cause of the age correlation here. More detailed chemical evolution 
modelling is clearly required to resolve the enduring puzzle of Ca 
abundances in passive galaxies. 
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Figure 9. The abundance ratios on the Fe-cr plane. The points show the Fe/H-cr relation for the combined Coma and Shapley galaxy sample, as in the left 
panel of Figure[2] The colour shading and contours were derived by heavily smoothing the X/Fe ratios onto this plane using a Kriging model (see e.g. Cressie 
1993). Contours are spaced at 0.05 dex intervals. Note the distinction between C and N, where the near-vertical contour indicate dependence only on ct, and 
Mg and Ca, where the contours are inclined and show a strong dependence on Fe/H. 



Table 2. The X-planes for the four elements Mg, Ca, C and N, derived from the combined sample of red sequence galaxies. The coefficients are derived from 
unweighted fits to the model [X/Fe] = ao + ai log (cr/100 km s -1 ) + 0,2 [Fe/H]. Bold face indicates coefficients with significance level >99%. For the 
orientation of the coefficient error ellipses, see Figure[8] For each ratio, we indicate the total dispersion around the fit, the contribution from the errors, and the 
estimated intrinsic contribution (all in dex). The scatter calculations include the effects of correlated Fe/H and X/Fe errors, and uncorrelated errors in a. 



ratio 


zero-point 


<r coeff 


Fe/H coeff 


total scatter 


error 


intrinsic 




ao 


ai 


a> 








Mg/Fe 


+0.12±0.01 


+0.35±0.03 


-0.33±0.04 


0.07 


0.06 


0.04 


Ca/Fe 


+0.04±0.01 


+0.15±0.05 


-0.40±0.05 


0.10 


0.08 


0.05 


C/Fe 


+0.06±0.01 


+0.43±0.04 


-0.10±0.04 


0.09 


0.08 


0.06 


N/Fe 


+0.10±0.02 


+0.41±0.05 


-0.07±0.05 


0.11 


0.10 


0.04 



Table 3. Extended multi-parametric fits to the abundance ratios, showing the generally non-significant coefficients of age and the offset term. The fits are to 
a model [X/Fe] = b + b\ log (cr/100 kms" 1 ) + b 2 [Fe/H] + b 3 log (t SSP /10 Gyr) + 6 4 5. Here 5 = for Coma galaxies and S = 1 for Shapley 
galaxies; its coefficient 64 is included to test for any constant offset between the two samples. Bold face indicates coefficients with significance level >99%. 
The total rms scatter, given in dex, is similar to that of the X-planes in all cases. 



ratio 


zero-point 


<r coeff 


Fe/H coeff 


Age coeff 


offset 


rms 




bo 


61 


62 


b 3 


64 




Mg/Fe 


+0.07±0.02 


+0.31±0.04 


-0.37±0.05 


-0.02±0.05 


+0.05±0.02 


0.07 


Ca/Fe 


-0.05±0.03 


+0.22±0.06 


-0.54±0.06 


-0.23±0.07 


+0.06±0.03 


0.10 


C/Fe 


+0.07±0.03 


+0.43±0.06 


-0.09±0.06 


+0.01±0.07 


-0.01±0.03 


0.09 


N/Fe 


+0.11 ±0.04 


+0.40±0.07 


-0.06±0.07 


+0.01 ±0.08 


+0.00±0.03 


0.11 
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Figure 10. Comparing abundance ratio data for milky way stars (points) and passive galaxies (lines and shading). The filled and open circles indicate thick 
and thin disk stars, respectively, from Bensby et al. (2005). In the N/Fe panel, squares show stars from Ecuvillon et al. (2004) and Israelian et al. (2004). 
In each panel, the X-plane fits (Table [2) have been used to derive the mean X/Fe-Fe/H relation for galaxies of velocity dispersion a = 50kms — 1 and 
a = 150 km s~ 1 (dashed and solid black lines, respectively). The agreement between the stellar data and the external galaxies suggests a degree of similarity 
in their enrichment processes. 



5 CONCLUSIONS 

We have analysed the correlations followed by the abundance ratios 
Mg/Fe, Ca/Fe, C/Fe and N/Fe, in a sample of red-sequence galax- 
ies over a wide range in mass and metallicity, using data of high 
signal-to-noise and uniform treatment. Using this sample we have 
discovered new systematic trends in the abundances. In particular, 
we have demonstrated that in general the abundance ratios X/Fe are 
driven by at least two partially decoupled factors, as evidenced by 
their planar relationships with Fe abundance and velocity disper- 
sion. 

We have shown that the a-abundance ratios Mg/Fe and Ca/Fe 
simultaneously decrease with Fe/H and increase with a. We inter- 
preted Fe/H dependence as due to a variation in the duration star- 
formation in different galaxies at given a, and hence different de- 
grees of Fe enrichment from SNIqS For both C/Fe and N/Fe, no 
correlation with Fe/H is observed at fixed a. This absence of cor- 
relation is consistent with these elements being produced primarily 
by low/intermediate mass stars, and hence on a similar time-scale 
to the Fe enrichment. We find a striking quantitative similarity in 
the X/Fe vs Fe/H trends (at fixed a), with the equivalent trends 



6 During final revision of our paper, Graves, Faber & Schiavon (2008) have 
confirmed the Mg/Fe-vs-Fe/H anti-correlation at fixed a in a sample derived 
from the Sloan Digital Sky Survey. 



for galactic stars, which suggests a high degree of regularity in the 
chemical enrichment history of galaxies. 

After allowing for the Fe/H effect, the Mg/Fe, C/Fe and N/Fe 
ratios all increase with a, with a similar slope. Coupled with the 
very different production sites and time-scales for Mg versus C and 
N, this is a powerful constraint on the mechanism responsible for 
the a correlations. Several traditional explanations, such as a sys- 
tematic variation of star-formation time-scale with a (e.g. Thomas 
et al. 2005), seem to be disfavoured at face value, because C and 
N are produced by low/intermediate-mass stars, and so should have 
time-scales similar to Fe. The shallow a dependence of Ca/Fe, in 
comparison to Mg/Fe, remains anomalous, as does the significant 
additional dependence on galaxy age, which is found only for Ca. 

In this paper, we have introduced a new diagnostic tool for de- 
scribing abundance ratios in galaxy samples, and made an initial 
exploration of the results obtained for four light elements. Our in- 
terpretation of the trends has been necessarily simplistic here, but 
the discussion indicates the greater power of the X-planes, in com- 
parison to the traditional one-parameter scalings typically analysed. 
We hope that the results will stimulate comparisons with more so- 
phisticated chemical evolution models, to understand better how 
the X-planes reflect the star-formation and enrichment history of 
red-sequence galaxies. 
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APPENDIX A: SEPARATE COMA AND SHAPLEY FITS 

In this appendix, we provide results for the single-parameter and 
X-plane fits obtained for the Shapley sample and the Coma sample 
taken separately. 

Table IAT1 gives the coefficients of the single-parameter fits. 
The "Shapey only" and "Coma only" columns correspond to red 
and blue fit lines in Figures [TJ [3] [5] and [6] In Table I A2] we present 
the X-pane coefficients for the two cluster samples separately, for 
comparison to the combined sample. The confidence intervals on 
these coefficients are compared in FigurefATJ 

These tests demonstrate that the X-plane signal is readily re- 
covered using the Shapley data alone, although the coefficients er- 
rors are larger, by up to a factor of two. This result further con- 
firms that the correlations reported for the combined sample are 
not driven by unidentified offsets between the two samples. Taking 
only the Coma galaxies, the errors are substantially larger, espe- 
cially on the a coefficients, because the range in a spanned by the 
Coma dwarfs is very narrow. The range of Fe/H however is fairly 
large, and the Coma sample alone is enough to recover the distinc- 
tion between Mg and Ca, which show significant anti-correlation 
with Fe/H, and C and N which do not. 
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Table Al. The single-parameter fits for X/Fe, for the combined sample, and for the two clusters fit separately. 







Combined sample 


Shapley only 


Coma only 


Ordinate 


X 


zero-point 


slope 


zero-point 


slope 


zero-point 


slope 


1 I a \ 
MlOOkms-O 


Mg 


+0. 17±0.01 


+0.18±0.04 


+0.l7±0.0l 


+0.14±0.05 


+0.21 ±0.04 


+0.2s±0.l2 




Ca 


+0.10±0.02 


-0.05±0.05 


+0.10±0.02 


-0.07 ±0.07 


+0.17±0.06 


+0.l2±0.16 




C 


+0.07±0.01 


+0.38±0.04 


+0.07±0.01 


+0.44±0.06 


+0.05±0.03 


+0.30±0.08 




N 


+0.11±0.01 


+0.37±0.04 


+0.ll±0.01 


+0.37±0.05 


+0.13±0.05 


+0.42±0.14 


( L \ 


Mg 


+0.11±0.01 


+0.05±0.02 


+0.l5±0.04 


+0.01 ±0.03 


+0.10±0.01 


+0.04±0.05 




Ca 


+0.12±0.01 


-0.03±0.02 


+0.l8±0.04 


-0.08±0.04 


+0.12±0.02 


-0.0l±0.06 




C 


-0.07±0.01 


+0.13±0.02 


-0.2 1 ±0.04 


+0.27±0.04 


-0.06±0.01 


+0.l3±0.04 




N 


-0.01i0.01 


+0.12±0.02 


-0.10±0.04 


+0.20±0.04 


-0.01±0.02 


+0.09±0.06 




Mg 


+0.18±0.01 


+0.28±0.04 


+0.20±0.0l 


+0.30±0.05 


+0.16±0.02 


+0.24±0.05 


Ca 


+0.14±0.02 


+0.1 3 ±0.05 


+0. II ±0.02 


+0.10±0.08 


+0.17±0.02 


+0.l9±0.07 




C 


+0.03±0.01 


+0.21±0.05 


+0.08±0.02 


+0.25±0.10 


-0.03±0.01 


+0.12±0.04 




N 


+0.07±0.02 


+0.19±0.05 


+0.l3±0.02 


+0.26±0.08 


+0.01 ±0.02 


+0.08±0.07 


[Fe/H] 


Mg 


+0.09±0.01 


-0.16±0.04 


+0.13±0.0l 


-0.24±0.06 


+0.00±0.02 


-0.28±0.05 




Ca 


+0.03±0.01 


-0.29±0.04 


+0.05±0.02 


-0.34±0.07 


-0.00±0.03 


-0.34±0.06 




C 


+0.02±0.02 


+0.11 ±0.04 


+0.06±0.02 


+0.12±0.10 


-0.07±0.02 


-0.04±0.05 




N 


+0.07±0.02 


+0.11 ±0.05 


+0.l2±0.02 


+0.17±0.08 


-0.03±0.03 


-0.06±0.07 



Table A2. The X-plane coefficients for the combined sample, and for the Shapley and Coma samples separately. Bold face indicates coefficients with signifi- 
cance level >99%. 





Combined sample 


Shapley only 


Coma only 


ratio 


(7 coeff 


Fe/H coeff 


a coeff 


Fe/H coeff 


a coeff 


Fe/H coeff 




ai 




«i 




«i 


a 2 


Mg/Fe 


+0.35±0.03 


-0.33±0.04 


+0.27±0.04 


-0.40±0.05 


+0.41±0.09 


-0.33±0.05 


Ca/Fe 


+0.l5±0.05 


-0.40±0.05 


+0.07±0.06 


-0.40±0.07 


+0.29±0.12 


-0.42±0.07 


C/Fe 


+0.43±0.04 


-0.10±0.04 


+0.48±0.06 


-0.12±0.08 


+0.33±0.09 


-0.07±0.05 


N/Fe 


+0.41±0.05 


-0.07±0.05 


+0.37±0.06 


-0.02±0.07 


+0.46±0.14 


-0.10±0.08 



Combined 



Shapley only 



Coma only 




0.2 0.4 
coefficient of log a 



0.2 0.4 
coefficient of log a 



0.2 0.4 
coefficient of log o 



Figure Al. Confidence intervals on the X-plane coefficients, for the combined sample and for the two clusters separately. Ellipses show the 68% and 95% 
confidence intervals on the two slope parameters. 
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